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^ : Abstract 

We calculate the range corrections to S-wave neutron-deuteron scattering in 
the doublet channel (5 = 1/2) to first order in r/a where a is the scattering 
, length and r the effective range. Ultraviolet divergences appearing at this 

order can be absorbed into a redefinition of the leading order three-body 
force. The corrections to the elastic scattering amplitude below the deuteron 
^ ■ breakup threshold are computed. Inclusion of the range corrections gives good 

■ agreement with measured scattering data and potential model calculations. 
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There has been much interest recently in applying Effective Field Theory (EFT) methods 
to nuclear physics [0-^- EFT provides a framework in which to exploit the separation of 
scales in physical systems in order to perform systematic, model-independent calculations. 
For nuclear few-body systems, the long-distance scale is set by the large two-body scattering 
lengths, while the short- distance scale is set by the range of the nuclear force. The EFT 
includes long-distance physics explicitly, while corrections from short- distance physics are 
calculated perturbatively in an expansion in the ratio of these two scales. In the two-body 
system for small momenta [p < m^r), this program has been very successful and calculations 
of deuteron properties and electroweak processes have been carried out to 1% accuracy (for 
a recent review, see Ref. 0]). 

In the nuclear three-body system, considerable progress has been made as well. In 
most three-body channels, the two-body EFT can be extended in a straightforward way 
However, the S-wave in the doublet channel (S* = 1/2) of neutron-deuteron scatter- 
ing (with the triton as a three-body bound state) is more complicated and exhibits some 
surprising phenomena [pHTTH. The renormalization of the three-body equations requires a 
one-parameter three-body force at leading order whose renormalization group evolution is 
governed by a limit cycle [|12|. The variation of the three-body force parameter gives a com- 



pelling explanation of the Phillips line (an essentially equivalent explanation was previously 
given in Refs. |1^,1|]). The phase shifts for S-wave neutron-deuteron scattering in the dou- 



blet channel have been studied at order (r/a)°, where a is the scattering length and r the 



effective range |T^. In this paper, we calculate the linear corrections in r/a to the elastic 
scattering phase shifts below the deuteron breakup threshold. The linear range corrections 
to the Phillips line have been studied in Ref. |T^ using a different formalism. 

Elastic neutron-deuteron scattering below the threshold of deuteron breakup can be 
described by an effective Lagrangian that includes only nucleons and has no explicit pions. 
For three-body calculations, it is convenient to use the dibaryon formalism of Ref. JTSf in 
which auxiliary fields with baryon number two are introduced to represent two-nucleon states 



in a given partial wave. The effective Lagrangian for the nucleon-deuteron system is [12 



C = NUtdo + ^)n- tUtdo + ^-AAti- slitdo + - aAs^ (1) 



2MJ ' V 4M 7 ""V 4M 
-^(tlN^T2aia2N + h.c.) - ^{slN^cr^r^r^N + h.c, 

-G^N^{g1{ti(yi)hvav + \gt9s \{ti(yi)^ SmTm + h-cl + fi'^(s„^r„)"^s^/r„/ + . . . 
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where N represents the nucleon field and ti (s^) are the dibaryon fields for the '^5'i {}Sq) 
channels and carry spin (isospin) one, respectively. The dots indicate higher order terms 
with more fields/derivatives, which do not contribute to the order we are working. The first 
line in Eq. (||) contains the kinetic terms for the fields A^, and Sm- The second line gives 
the coupling of the dibaryon fields to nucleon fields where ai (r^) are Pauli matrices acting 
in spin (isospin) space. Finally, the third line contains the three-body force. Note that this is 
the only three-body operator without derivatives that preserves spin and isospin symmetry 
1^,16]. Other non-derivative three-body operators can be related to the one shown via Fierz 



transformations. Since the action is quadratic in the fields ti and s^, it is straightforward to 
integrate them out and show that the theory is equivalent to one of nonrelativistic nucleons 
interacting via two-body and three-body contact interactions |]7|,|8|JT^. 
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FIG. 1. Geometric series leading to the exact dibaryon propagator. Double (single) lines rep- 
resent bare dibaryon (nucleon) propagators, respectively. 

To obtain the exact dibaryon propagator, the bare propagator must be dressed by nucleon 
loops to all orders. This is illustrated in Fig. |1]. The diagrams form a geometric series and 
are easily summed; the result is 

iDj{qo, q ) = — , , (2) 

-7j + (7I + Mqo - q^/A) + y'-Mgo + g 74 - 

where the subscript j = t (s) for the ^Si{^So) channel. This propagator has a pole at 
Mqo — = — 7|, where 74 = 45.68 MeV and 7^ = —7.88 MeV. The effective ranges 
are rt = 1.76 fm and = 2.75 fm. These effective range parameters are related to the 
parameters appearing in the Lagrangian via 



7i = 



Wn^Ai \ , ATT 



The scattering length is given by —1/aj = — 7^ + rj7|/2. Thus the Lagrangian of Eq. (|1]) 
reproduces the first two terms in the effective range expansion of the two-body scatter- 
ing amplitude. Note that including the kinetic terms for the dibaryon fields in Eq. (|1|) is 
necessary to obtain the range correction. 

The power counting in an EFT makes it possible to organize calculations in a system- 
atic expansion in a small parameter. In the two-body sector, the power-counting scheme 
of [0,0 takes p ~ 1/a ~ Q, where p is the typical momentum of a nucleon and 1/a is 



the inverse scattering length. Note that 7 is 0{Q) as well since 7 = 1/a + 0{r/a?'). The 
expansion parameter of the theory is Q/A, where A represents the scale where short- distance 
physics becomes important. For systems interacting via short-range interactions, the effec- 
tive range is expected to be set by A. So the expansion parameter is Q/A ~ r/a ^ 7r.[| 
The leading order term in the expansion of the two-body scattering amplitude is 0{Q~^). If 
we integrate out the dibaryon fields, the two-body operators without derivatives are treated 
nonperturbatively because the renormalization group equations dictate that their coefficients 
are 0{Q~^). In the dibaryon formalism, 27rA/(M(7^) = 1/a ~ Q, which requires summing 
the bubbles in Fig. |T]. At leading order in Q, the theory reproduces Eq. (^) in the limit 
TtiTs — > 0. The two-body operators with two spatial derivatives or one time derivative 
which are given by the kinetic terms of the dibaryon fields appear at next-to-leading order 
in the Q expansion, giving a contribution of 0{Q^) to the scattering amplitude. When 
these operators are included in the three-body problem, they also give a correction that is 
suppressed by one power of Q relative to the leading order. Since the coefficients of these 



^Note, that it is sufficient to take ''ytft = Isfs = If for power counting purposes. 
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FIG. 2. Leading order integral equation. Shaded (full) double lines indicate spin-triplet 
(spin-singlet) dibaryon. Single line indicates nucleon propagator. 



operators are proportional to the two-body effective ranges, we refer to these corrections as 
effective range corrections. 

Next we consider power counting in the three-body sector. EFT power counting shows 
that all diagrams that contain only non-derivative contact interactions are 0{Q~^) ||12|| . 
These diagrams can be summed using the integral equation shown in Fig. |^. At this order 
in the EFT, it is appropriate to use the propagator of Eq. (^ in the limit rf,r<j 0. 
However, it turns out that the resulting integral equation has no unique solution if the 
cutoff A is taken to infinity |TB|JT^. If the integral equation is regularized with a finite 
cutoff A, the solution displays a strong cutoff dependence. The integral equation can be 
renormalized by absorbing the cutoff dependence into the three-body force of Eq. (|ip |]T2 . 
The diagrams with the three-body force are naively 0{Q^). However, in the 5 = 1/2 channel, 
the leading three-body operator with no derivatives is relevant at low energies because the 
renormalization group evolution enhances the coefficient G3 in Eq. (|l]) to 0{Q~^) rather 
than 0(Q°) as naively expected ||12|.0 It is convenient to pull out a factor of 2M/A^ and 
define G3(A) = 2MH{A)/A^ with 



HiA) 



sin[so ln(A/A^) - arctg(l/so)] 
sin[soln(A/A*) + arctg(l/so)] ' 



(4) 



where Sq ~ 1.0064 is determined by the asymptotic behavior of the integral equation and 
A* is the three-body force parameter |]12|,^. As a consequence, there are certain cutoffs for 
which the three-body force vanishes. Since all observables are independent of the cutoff, it 
is possible to obtain a renormalized equation by choosing a cutoff with vanishing three-body 
force 1^ . The parameter A* then appears in the upper limit of the integral. Evaluating the 
diagrams in Fig. |], the renormalized equation in the limit r^, — > takes the form 



12,19 



-7t + ppV^ - ME, 

p2 _ ^2 



aUp) = Kk{p, k) + - 

TT Jo 



2 M^A,) g2p 



dq- 



-7s + a/3pV4 - MEk , 
2 -blip) = SK.ip, k) 



— k"^ 

2 /-AnCA.) g2p 



Kk{p.q)(al{q) + ?>hl{q)) (5) 



TT JO 



' q^ — k"^ 



Kk{p,q)[3aUq) + bl{q) 



^Note that the three-body force in Eq. (|l]) does not contribute in the quartet channel because of 
the Pauli exclusion principle. 
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where k (p) denote the incoming (outgoing) momenta in the center-of-mass frame, MEk = 
3/c^/4 — is the total energy, and A„ = A,,exp[(n7r + arctan(l/so))/so] with n a natural 
number. Three-body observables are independent of n up to corrections that are suppressed 
by inverse powers of A„ The kernel Kk{p,q) arises from the S-wave projected one- 



nucleon exchange and is given by 



Kk{,P, q) = —\Ti\ —- — — . 6) 

2pq \q'^ — pq + p'^ — MEk 



The amplitude a°(p) is normalized such that a^lk) = l/(A;cot 6) with 5 the elastic scattering 
phase shift. 

Recently, various authors have suggested treating range corrections nonperturbatively 



in both the two- and three-body systems [22-26]. This can be motivated by arguing that 



a nonperturbative treatment of range corrections resums large corrections proportional to 
iltft)"' to all orders in n. Since jtft ~ 0.4, these corrections can be numerically important 
despite being formally subleading in the Q expansion ||2^. Alternatively, one can imagine a 
power counting in which r is 0{Q~^) [|^-^. In the two-body sector, such a power count- 



ing was shown to simplify calculations and improve the convergence of the expansion p5 
Furthermore, the coefficients of higher derivative S-wave operators are no longer enhanced 
by renormalization group evolution and naive dimensional analysis can be used to estimate 
their contribution to amplitudes. 

A nonperturbative treatment of range corrections in the three-body problem was studied 
in Refs. [|^,^. The integral equation of Fig. || is solved using the propagator of Eq. (H) 
without expanding in r. This drastically changes the nature of the solution to the integral 
equation. Since the dibaryon propagator falls as rather than 1/q for large q, the kernel 
is damped at large loop momenta, and the integral equation has a unique solution even in the 
absence of the three-body force. The three-body force is not enhanced by renormalization 
and is a subleading effect suppressed by Q^. There is no three-body parameter in the 
leading order calculation and therefore no Phillips line. Surprisingly, the Phillips line is 
not even recovered at 0{Q^), when the three-body force is included. In Ref. [^, it was 
shown that when the range is treated nonperturbatively and the cutoff. A, is taken to 
infinity, the solution is completely insensitive to the numerical value of the three-body force. 
Furthermore, the obtained scattering phase shifts strongly disagree with experiment PB| . 

It is not clear why the nonperturbative treatment of the range corrections fails. One 
possible problem is that the dibaryon propagators have spurious poles at Mq^ — q'^/A = 2/r + 
0(7r). These poles can be avoided by using a momentum cutoff in the range 7 < A < 1/r. 
The three-body force can be tuned in such a way as to leave results approximately cutoff 
independent when the cutoff is varied within this window. Unfortunately, since at ~ 3rt, 
unreasonably small cutoffs must be used in this method. Moreover, the nonperturbative 
range correction worsens the agreement with the Phillips line. 

In the present paper, we take a different approach and compute the range corrections to 
a^(A;) perturbatively, working to 0{Q~^). An important question is whether higher derivative 
three-body operators also contribute at this order. These operators will include at least one 
time derivative or two spatial derivatives and hence contribute to the amplitude at O(Q^), 
according to naive dimensional analysis. This is much higher order than the range correction 
to the three-body amplitude, which is 0{Q~^). Higher derivative three-body operators 
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FIG. 3. Feynman diagrams for range correction, (a) Range correction: double line with cross 
denotes O^jr) piece of the dibaryon propagator, (b) Subleading three-body force: filled square 
denotes insertion of H^{A). The propagators are as in Fig. |2|. Not shown are diagrams that vanish 
as A ^ oo. 



could only contribute at this order if there were a renormalization group enhancement of 
their coefficients, which is not the case. Consequently, the range correction is the only 
contribution at this order. In the following, it is more convenient to label the contributions 
by powers of 7r relative to the leading order. The linear range correction is then 0(7r). 
Below, it is demonstrated that the range correction can be calculated without introducing 
any terms not already present in Eq. (|1|). The operator which renormalizes the leading order 
calculation can also be used to renormalize the loop graphs that appear at 0(7r). However, 
the running of the coefficient H{A) needs to be modified. 

Since the computation of the range corrections requires no additional counterterms, only 
the two-body effective ranges and enter as new parameters at this order of the calcu- 
lation. The only parameter to be fixed from a three-body datum is the three-body force 
parameter, A^,. For example, if A^, is fixed to reproduce the observed neutron-deuteron scat- 
tering length, the energy dependence of the scattering amplitude is completely predicted up 
to corrections of 0((7r)^). 

To calculate the range corrections it is convenient to formally expand the amplitudes 
0'k{p)i bk{p) and the three-body force H{A) in powers of •yr, 



b{q)=bl{q)+bl{q) + ... 
H{A) = H\A) + H\A) + ... 



(7) 



where the superscript denotes the power of 7r. The amplitudes a1{q) and b1{q) are the 
solutions of Eq. (|^) and H^{A) is given in Eq. (^). In the following, we will calculate al{k) 
and obtain the renormalization group evolution of H^{A). 

The Feynman diagrams contributing to al{k) are shown in Fig. ^ Fig. ^(a) shows the 
diagrams from the range correction. The double lines with the cross denote the 0(7r) piece 
of the dibaryon propagator, Dj{qo, q). Expanding Eq. (0) in r^, we obtain the Feynman rule 
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^DA%^ ? ) = T7-r 77 ^ , ' =-2 (8) 



In order to evaluate the Feynman diagrams in Fig. ^(a), it is convenient to use the form given 
in the second hne of Eq. (j^). The constant term in Eq. (|^) does not contribute because the 
contour integral over the zero component of the loop momentum can be evaluated without 
enclosing a singularity |2^. Fig. |^(b) gives the contribution of the subleading piece of the 



three-body force, H^{A). In addition to the diagrams shown in Fig. ^(b), there are four more 
diagrams where the -ff^(A) insertion is dressed on only one side by either aP^^q) or b\{q) and 
a diagram with the undressed H^{A). These diagrams are also included in the calculation, 
but it is shown below that their contribution vanishes as A oo. 

The evaluation of these diagrams is straightforward. After projecting onto the S-waves, 
we obtain: 

W) =-itrtat^[k)+-it- dq— — (9) 

^ -^0 9^ - -it + ^J?>q^/4.- MEk 

/„ "^w^' r-'-- h + Av?^ 

2 M„(A.) q2p 



3 AS 



1 + -/ dq^^ial{q) + hl{q) 

Ti Jo q"^ — k"^ \ 



Note that for the integral in the second line no principal value prescription is required since 

m) = 0. 

There are four contributions to the 0{'-fr) corrections to the amplitude. First, the residue 
of the pole in the deuteron propagator changes by a factor 1 + jt'^f This changes the LSZ 
factor that goes into the leading-order calculation, giving rise to the first term in Eq. @. 
The second and third terms in Eq. (H) come from the diagrams in Fig. 1(a). The fourth 
term comes from diagrams involving the subleading piece of the three-body force, H^{A). 

In order to determine H^{A), it is necessary to know the large A dependence of the 
integrals in Eq. (|^). Since the combination a°(g) — bl{q) falls off sufficiently fast for large 
q (cf. Ref. ||12[|), only the asymptotic form of the combination a^(g) + bl{q) for g 7 is 



needed: 



+ bliq) ^ Skiq) = U{k) cos{s, ln(g/A,) + 5) , (10) 



where we have suppressed the dependence of Nik) on A and 7.^ It is especially important 
to note that the phase of Sk{q) is independent of k so Sk{q) factorizes into a product of 
an unknown function of k and a known function of q. This will allow us to determine 



^This is can be done safely because N{k) is of 0(A°) H^. 
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iJ^(A) analytically and show that all divergences in the range correction are cancelled by 
this counterterm. 

Inserting Eq. (^) into Eq. (j^) gives the divergent piece of the range correction 



_ 27t / Itn + IsTs dq 

oa^iv — 7=- 

TT \ V3 J q 



(11) 



The divergence is cancelled if 



TT 



A2 



16 
1^ 



(7tn + IsTs) 



/dqsl{q)/q 
[S''dqs,{q)Y 



(12) 



+ 25 + sin(2so ln(A/A,) + 25) 



4soA2 

F(A) 
A2 ' 



[cos(soln(A/A*) + 6) + sq sin(so ln(A/A*) + 6)Y 



where the last line defines F{h). Note that -ff^(A) is the coefficient of an operator with 
no derivatives and cannot be a function of k. As required, all k dependence has cancelled 
in the expression for if^(A). Furthermore, H^{h)/h? scales like (up to logarithmic 

corrections) for large A. In the third line of Eq. (|^), it is only necessary to keep the term in 
which H^{A) is multiplied by two linearly divergent integrals. This term corresponds to the 
diagrams shown in Fig. ^(b). The other terms can be discarded since their contribution can 
be made arbitrarily small by choosing an appropriate value for A. Thus the renormalized 
expression for the range correction is 



a.ik) ^ 7.n a,ik) + 7.- dq^ ^ ^3.2/4 - ME. 

2 .A.(AO qHljqy 8 
^^Vio ^(g2-A;2)2 3^^^^ 



(13) 



27t . ^ .n^n) 



TT 



n 



A„(A,) 



-7. + V3g2/4 - MEk 
\ 1 2 

dq(al{q) + hl{q) 



The above expression is cutoff independent (up to corrections of 0(1/A„)). Note the func- 
tion -F(A) is known analytically but depends on the asymptotic phase 5, which must be 
determined numerically by fitting to the leading-order solution a^(g) + b^{q). 

In Fig. ^ we show our result for k cot 6. The leading-order (LO) calculation is indicated 
by the full line, while the next-to-leading-order (NLO) calculation is given by the dashed 
line. At each order, A* is tuned to produce the measured neutron-deuteron scattering 

We find Af'^ = 3.6 7f = 0.83 fm"^ and A^^^ = 



length, a^^P 



(0.65 ± 0.04) fm |28| 



4.1 7t = 0.95 fm . In Fig. ^, the circles correspond to the phase shift analysis of Ref. p9 



while the squares show a potential model calculation using the Argonne V18 nucleon-nucleon 



potential and the Urbana three-nucleon force |30]. The triangle gives the experimental value 
of the nd-scattering length from Ref. |2^. The range corrections are small all the way up to 
the breakup threshold. It is encouraging to see that the perturbative corrections are small 
even though 'jtTt ^ 0.4 is not a very small expansion parameter. This suggests that the 
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FIG. 4. Results for kcot6: LO (NLO) calculation is indicated by solid (dashed) line. Triangle 
gives the experimental value of the scattering length from Ref. |28|, circles show the phase shift 
analysis of Ref. , and squares show the potential model calculation of Ref. . 



EFT expansion is well behaved. The range correction clearly improves agreement with the 
phase shift analysis of Ref. [^]. Note that this phase shift analysis is more than 30 years 
old and gives no error estimates. The errors of the analysis are at least as large as the error 
of the scattering length; most likely they are larger. Consequently, it is more meaningful to 
compare to the potential model calculation of ||3^, which agrees well with the NLO result. 

In summary, we have calculated the S-wave phase shifts for neutron-deuteron scattering 
in the doublet channel to 0{'~fr) and found good agreement with available data. We have 
shown that the corrections at this order can be renormalized by modifying the running of the 
leading order three-body force. Apart from the two-body effective ranges and r^, no new 
parameters enter at this order. In Refs. it was shown that the perturbative treatment 
of the range corrections in other channels where three-body forces are subleading gives good 



agreement with available data as well. As stated earlier, Ref. [0] has calculated the Phillips 
line to 0(7r) and obtained results which agree well with the Phillips line obtained from 
various potential models. Furthermore, if one demands that the neutron-deuteron scattering 
length be correctly reproduced, then the Phillips line predicts the triton binding energy. The 
prediction is 8.0 MeV at 0((7r)°) |12|,|1|1 and 8.8 MeV at OM 0. These numbers compare 
well with the measured binding energy of 8.5 MeV, and again the range corrections improve 
agreement with experiment. Together, these results show that the power counting of |T7,TE 
is adequate for three-body nuclear systems at very low energies. This suggests that the 
perturbative method could be used for precise calculations of phenomenologically important 
three-body processes such as polarization observables in neutron-deuteron scattering and 
the /?-decay of the triton. 
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